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We propose a scheme for ultrafast real-time imaging of laser-induced collective electron oscillations (Mie plasmons) in gas phase 
metal clusters by interferometrically stable scanning of two intense few-cycle optical laser pulses. The feasibility of our nonlinear 
spectral interferometry method with experimentally accessible observables is tested in a theoretical case study on simple-metal 
clusters (Nai47). The results show that the plasmon period and lifetime as well as the phase and relative amplitude of the collective 
ry j electron motion can be extracted with sub-fs resolution. The access to nonlinear response effects, as the demonstrated increase of 
the plasmon lifetime with laser intensity due to ionization-induced contraction of the electron cloud, opens up vast opportunities 
f° r interrogating ultrafast many-particle dynamics in nanosystems under strong laser fields with unprecedented resolution. 
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1 Introduction 

Coherent control of ultrafast processes in atoms and molecules 
has become a vital research branch of present-day laser-matter 
science!^. Nowadays, the availability of phase-controlled 
few-cycle pulses and shaped light fields with precisely ad- 
justable amplitude, phase, and polarization^"^ offers a ver- 
satile toolbox for probing and steering electronic dynamics 
down to sub-fs time scales. Such tailored light fields have 
been applied for illuminating various aspects of ultrafast non- 
linear laser-matter interactions such as quantum path control 
for generating selected high harmonics* 11 ! 12 ! or single attosec- 
ond pulses 13 14 , tunneling of electrons in atomic photoioniza- 



ti on LL5J va i ence electron motion in atoms^H electron localiza- 
tion in molecular dissociation 17 , and selective population of 
dressed electronic states^. 

So far, control schemes with sub-fs resolution have been 
applied predominantly to processes, whose main physics can 
be understood in an effective single active electron picture. 
A promising and complementary future direction of ultrafast 
laser science is the manipulation and imaging of many -particle 
electron dynamics in the nonlinear response domain. From 
such efforts, deeper insights into the nature of collective elec- 
tronic processes in ultrashort laser fields 19 " 21 ! their ultrafast 
control 22 , and potential new applications can be expected^!. 

An instructive example for such many-particle dynamics 
are collective electron oscillations (Mie plasmons) in simple 
metal clusters^. In a simple picture, the plasmon can be con- 
sidered as a coherent oscillation of the valence electron cloud 
against the charged ionic background in a spherical metal- 
lic drop. Within the limit r c i uster <C A, only the dipole mode 
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couples to the laser field and may be excited resonantly with 
high cross sectionEKz] a simple estimate for the spectral 
position of the Mie plasmon 28 for a spherical metallic drop 
is 0>Mie = \/e pi/3£om e , which solely depends on the charge 
density of the ionic background pu while £n, e, and m e denote 
the vacuum permittivity, the elementary charge and the elec- 
tron mass. In reality, additional effects like electronic spill 
out, polarization of core electrons, electronic quantum con- 
finement, and cluster deformation lead to shifts and spectral 
fragmentation of the collective resonance^E] For our pur- 
poses, however, the plasmon can be considered as a single 
resonance with an energy ft (Ores of a few eV, typically well 
above the photon energy of near-infrared (NIR) laser fields, 
and a lifetime T res of a few fs. 

In the weak field regime, the spectral position, width, and 
oscillator strength of the plasmon in gas phase clusters may be 
determined by spectral measurements of photofragmentation 
or ph otoioniza tion cross sections with tunable narrow-band 
lasers 24 l 26 l 27 l 32 l For deposited cluster also two-photon pho- 
toemission has been applied^. When aiming at the strong- 
field domain, a time-domain analysis of plasmon oscillations 
with sub-cycle resolution becomes desirable, as ionization and 
cluster explosion result in transient optical properties. This 
has not been achieved yet, though the time-dependence of the 
plasmon mode itself is routinely exploited in strong-field ex- 
periments. For example, the expansion-induced red-shift of 
the plasmon frequency occurring on the time scale of nuclear 
motion makes it possible to establish a plasmon reso nance 
with NIR laser pulses, e.g. in a pump-probe setup 3EM2D a 
sub-fs analysis of plasmon oscillations may thus be a chal- 
lenging stepping stone towards a deeper understanding of col- 
lective processes in intense laser fields. 

In the present study we propose a possible route towards 



This journal is ©The Royal Society of Chemistry [year] 



Journal Name, 201 0, [vol], 1 -[8] | 1 



such ultrafast time-resolved analysis. To this end we develop 
and theoretically explore a scheme to excite and image plas- 
mon oscillations in expanding simple metal clusters with sub- 
fs resolution by nonlinear spectral interferometry with few- 
cycle laser pulses. The excitation and analysis of the plasmons 
is utilized by a pair of few-cycle pulses, where the plasmon 
oscillations are induced by the first pulse and probed with the 
second one. This enables a time-resolved mapping of the co- 
herent electron motion by scanning the delay between the few- 
cycle pulses with interferometric accuracy. The performance 
and accuracy of this approach is demonstrated with semiclas- 
sical Vlasov simulations on the model system Nai47, where 
we show the possibility to resolve the dependence of the plas- 
mon lifetime on the intensity of the pump pulse, which is a 
truly nonlinear effect. A key finding of our analysis is that 
experimentally accessible observables like ionization or pulse 
depletion can be used for the analysis, even if the system re- 
sponse is nonlinear. 

The manuscript is structured as follows. Sec. [2] introduces 
the strategy and approximations of our spectral interferome- 
try approach and presents benchmark results for a simplified 
oscillator model. Sec.[3]discusses the semiclassical Vlasov ap- 
proach and presents the numerical results on the time-resolved 
imaging of plasmon oscillations for representative cases. Con- 
clusions and an outlook are given in Sec. [4] 



2 Spectral interferometry 

Spectral interferometry is the basis of Fourier transform spec- 
troscopy and has also been recognized as a powerful tool for 
coherent control of molecular reactions with light 37 . It has 
been successfully applied to various scenarios, such as the 
control of electronic wave packet motion in atomsEHE] rov i- 
brational nuclear dynamics^, and quantum interferences in 
molecules 41 . 

For a brief introduction of key aspects relevant for the ap- 
plication of spectral interferometry to collective processes and 
the derivation of important approximations, a simplified case 
in linear response is considered. A suitable benchmark sce- 
nario is the excitation of a damped harmonic oscillator by a 
pair of linearly polarized few-cycle laser pulses. Based on 
that, we discuss and test two approximations for the total en- 
ergy absorption as function of pulse delay, which are key to 
the retrieval of the to-be- analyzed dynamical properties. 



2.1 Basic equations 

We begin with the electric field of a single pulse 



E (t) = ]:E em (t)e 



where (Oh is the angular frequency of the carrier wave, (Pce 
is the carrier-envelope phase (CEP) and E Qm (t) is the real- 
valued temporal field envelope. The total field of the pulse 
pair E(t) = Eo(t) + Eo(t — At) then reads 

p —i(pCE r -i 

E(t) = [^(O^+^env^-AO^^] +C.C, 

(2) 

where At is the pulse delay. Note that we require pulses 
with the same CEP, which may be realized by a Michelson 
or Mach-Zehnder interferometer setup. The total field in the 
Fourier domain E((o) = J^E^e^dt has the form 



E(co) = 



1 



JcoAt 



[e- i( P^E em (co -gjq) + e i( P^E em (co + coq) 



(3) 

where E mv (co) is the Fourier representation of the envelope. 
The net energy absorption of the oscillator driven by the light 
field E(t) can be expressed in the Fourier domain as an integral 
over the product of the spectral intensity 1(g)) = 1/271 \E(g))\ 2 
and the absorption cross section of the oscillator <j((o) via 



W abs (At) 



i_ r 

2% J-c 



E(CQ)E*(CQ)o(CQ)dCQ. 



(4) 



Assuming E em (co — (0n)£* nv ((0 + (0 o) = 0, i.e. a sufficiently 
well peaked spectrum of the envelope function, and exploiting 
that a(co) = o(— co) for a classical oscillator we can write 



W abs (At) : 



j_ r 

271 J-c 



( 1 + cos coAt) \E em (co - coo)\ 2 a(co)d(0. 



(5) 

By this reformulation all spectral contributions around 
CO = — C0o have been shifted to positive frequencies. It should 
further be noted that the CEP has dropped out. For sufficiently 
weak damping of the oscillator we can now express the cross 
section by a normalized Lorentzian and obtain 

(1+cos coAt) \E Qm (co-coo)\ 2 /27i 



Wabs(Af) 



71 J-c 



-dco, (6) 



-z(o) ?+<Pce) . 



-c.c, 



(1) 



((Ores- CO) 2 + J 2 

where (Ores and y specify the frequency and spectral width (in- 
verse lifetime) of the oscillator resonance. The numerator in 
the integral is the intensity spectrum of the double-pulse laser 
field I(co). In this formulation of the spectral intensity all sig- 
nificant spectral contributions have been shifted to the range 
of positive frequencies. The usual relation f™ oo I(G))dG) = 
f™ OQ E 2 (t)dt is still valid. The double-pulse spectrum 1(g)) 
contains the intensity spectrum of a single pulse 

I (G)) = \E &m (G)-G)o)\ 2 /47l (7) 

times a spectral modulation with 2(1 + cos G)At). This tunable 
interference of the pulses in the Fourier domain (cf. insets of 
Fig. [T^) is the key feature exploited in spectral interferome- 
try. For extracting the deeper physical meaning, the integral 
in Eq. ^ will now be considered in detail for two approxima- 
tions. 
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2.2 Linear spectrum approximation (LSA) 

Assuming an oscillator resonance well inside the laser spec- 
trum (detuning up to the spectral width of the laser pulse), 
Iq(co) can be linearized around (Ores by 



7^ SA (co) = a + (a)-a) res )j3, 



(8) 



with a = /o(cHfe s ) and j8 = d/o(fi>) / 'd(o\co=^- Carrying out 
the integral in Eq. ([6]) within this linear spectrum approxima- 
tion (LSA) yields an absorption of 



W, 



LSA/ 



abs 



(At) = 2a 1 + J\ + c£ SA cos (^Ar + (Plsa) e yAt 



(9) 

where c LS a tan^sA = Py/a. 

In the special case of an oscillator resonance at the maxi- 
mum of Iq((d), i.e. for /3 = 0, this reduces further to 



which contains the single pulse fluence Fq = i J^E^^dt 
and the autocorrelation of the envelope function A env (Ar) = 
J^Eqw^Eqw^ — At)dt. The term in square brackets in 
Eq. ( p3| can be identified with the fluence F(Af) of the double 
pulse laser field. Hence, a non-resonant excitation results in 
an absorption signal proportional to the field fluence, which 
can be exploited for identification of such scenario. 

2.4 Simple examples with gaussian few-cycle pulses 

We now demonstrate the performance of the above ap- 
proximations for representative benchmark examples with a 
damped oscillator. Considering gaussian pulses, the envelope 
function can be written as E Qm (t) = Ee~ l /T , where E is the 
single pulse field amplitude and T is the pulse width. The lat- 
ter can be related to the full width at half maximum of the 
intensity via T = TFWHM/V21n2. With the envelope function 
in the Fourier domain E em (co) = ET\f%e~® T / 4 the spectral 



W™ A (At) = 2a (l + cos (ofesAf ) e~ rAt ) ■ (10) intensit y of the sin § le P ulse is 



Hence, the delay-dependent absorption exhibits a static off- 
set plus a damped periodic oscillation with the eigenfrequency 
and lifetime of the oscillator. Most importantly, for vanishing 
CEP the oscillations even provide a direct time-domain im- 
age of the dipole velocity of the oscillator after the first pulse, 
which turns out to be the key to our imaging scheme. Apart 
from a phase shift by (JOlsa and a scaled modulation amplitude 
(both can be compensated for a known laser spectrum), this 
picture even persists if the resonance is located in the wings 
of the laser spectrum (/3 ^ 0), cf. Eq. (|9j. It will be shown 
later that these modifications are small even for considerable 
detuning of the order of the spectral width of the laser pulse. 

2.3 Constant cross-section approximation (CCA) 

As the second limiting case we consider a fully off-resonant 
excitation, i.e. co^s far outside the laser spectrum. The cross 
section can then be approximated with the value at the carrier 
frequency Go = g(q)o). The absorption integral within this 
constant cross section approximation (CCA) becomes 



/o(fi>) : 



-((0-(Oo) 2 T 2 /2 



(14) 



which yields clsa = tan (JPlsa = — ^ 2 y(&>res — (Do) for the linear 
spectrum approximation in Eq. ([9]). 

In the next step, optical excitation with few-cycle pulses 
(^fwhm = 3fs) at 800 nm carrier-wavelength (h(0o=l .54 eV) 
with zero CEP is considered. Exemplarily, we study three 
different oscillator configurations with resonances (Ores = 
1.3a)o, and 2(0$. For all cases a fixed damping constant of 
hy= 0.13eV (this corresponds to a lifetime T res = 5fs) is as- 
sumed. These examples reflect (i) fully resonant, (ii) near- 
resonance, and (iii) nonresonant excitation. Normalized cross 
sections are shown Fig.[T^ in relation to the single pulse inten- 
sity spectrum. 

Lets now consider some specific observable A, which may 
be absorption, mean electron energy, or electron yield in a real 
cluster experiment. The delay dependence of A(At) can then 
be expressed by the dimensionless modulation parameter 



M A (At) 



A(At) 



W^ A (At)=2G [ (1+ cos co At)Io(co)dCQ. (11) 

J — oo 



A(At' > X res ) 



1, 



(15) 



Applying the convolution theorem one finds 

/oo 
El w {t)dt (12) 
-oo 

+ CTocosa^A^ / E Qm (t)E Qm (t — At)dt 

J — oo 



Go [2F + A em (At) cos COoAt] 

V v ' 

F(At) 



(13) 



which vanishes for pulse delays much larger than the oscillator 
lifetime. The modulation parameter for the absorbed energy 
with the full expression in Eq. ^ then follows as Mw(At). 
This function is used to benchmark the modulation parame- 
ters M^f A (At) and M^ CA (At) obtained from the approximated 
energy absorptions in Eq. ^ or Eq. ( [TT] ), respectively. The re- 
sulting evolutions of My/ {At) as function of pulse delay (solid 
red curves) are shown in Figs.[T]3-[T]i in relation to the applica- 
ble approximation M^f A (At) or M^ CA (At) (dashed, approxi- 
mation as indicated). 
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Fig. 1 Spectral interferometry results with gaussian 3 fs pulses 
(hcOQ = 1.54eV) for three Lorentz oscillator settings, (a) single 
pulse spectral intensity (gray) and normalised oscillator absorption 
cross sections for o^es = CQo, 1.3ft)o, and 2(Oq (red, green, blue) for 
7 _1 = 5fs lifetime; insets show electric field evolutions (black) and 
resulting spectral intensity distributions (red) for different delays; 
(b-d) delay-dependent modulation signals Mw(At) (red) compared 
to results in linear spectrum approximation (LSA) or constant cross 
section approximation (CCA); corresponding (p^SA values as 
indicated; dashed blue curves show the pump-only dipole velocity 
signal Vd(f); insets depict evolutions of the laser field (black), dipole 
velocity (blue), and absorbed energy (red) for indicated delays. 



For the fully resonant and the near-resonant scenario 
(Figs. [lj> and [it), the predictions of the linear spectrum ap- 
proximation M\f A (At) (rescaled in amplitude as indicated) 
are in good agreement with the exact results Mw (At) and show 
damped oscillations with high contrast. A closer comparison 
of the rescaled M^f A (At) signals with Mw(At) shows, that the 
modulations appear in phase and with the oscillator frequency 
and decay time for At > 6 fs, i.e., outside the region with strong 
pulse overlap. Deviations for smaller delays reflect signal dis- 
turbance by higher order interference terms that are negligible 
outside the pulse overlap region. 

For selected delays, time-domain evolutions of the laser 
field, dipole velocity, and absorbed energy are sketched as in- 
sets and provide a direct image of the dynamics that is mapped 
into the modulation parameters. Focussing on the upper time- 
domain plot in Fig.fTJ), a resonant dipole oscillation is excited 
by the pump pulse. The fully resonant nature is reflected by 
the fact, that the dipole velocity and the laser field are in phase. 
For the chosen delay, the probe pulse arrives in phase with the 
dipole velocity (constructive interference) and the oscillator 
motion can be strongly amplified. The energy gain from the 
probe clearly exceeds the absorption from the pump pulse due 
to advantageous coherent superposition of remaining coherent 
pump-induced motion and the probe excitation. 

In contrast to that, the lower cartoon in Fig.fTJ) shows an ex- 
ample for destructive interference, where the probe pulse even 
stops and reverses the oscillator motion and the net energy 
gain from the probe is substantially reduced. The oscillation 
of the modulation parameter with pulse delay thus reflects the 
alternation of amplification and suppression of energy absorp- 
tion from the probe pulse due to the residual oscillator motion. 
In turn, the decrease of the modulation amplitude with delay 
indicates the decay of the pump-pulse-induced oscillator mo- 
tion. For long delays (At ^> T res ) the interaction with the probe 
pulse is incoherent, as no phase- sensitive excitation remains at 
the time of its arrival. 

In the fully resonant scenario (j3 = 0) in Fig. [TJ), the mod- 
ulation signals directly image the pump-induced motion of 
the oscillator. This can be inferred from the close-to-perfect 
agreement of M w (At) and the rescaled M^f A (At) with the 
dipole velocity v&(t) for pump-only excitation (blue dashed 
curve in Fig. [TJ)). Outside the pulse overlap region, signals 
are in phase and show the same relative amplitude evolution. 
This direct mapping of oscillator motion into the modulation 
parameter is the heart of our imaging scheme. 

A similar set of time-domain examples is shown for the 
near-resonance scenario in Fig. [T]: (see insets). The general 
trends and the modulation effects are very similar to the fully 
resonant case. It should be emphasized that the modulation 
signals (Fig. [T£) are now phase shifted to the pump-induced 
dipole velocity by a small phase offset of the order of (JOlsa- 
The latter could be determined by fitting the LSA formula 
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from Eq. ^ to a given Mw(At) signal outside the pulse over- 
lap region or directly from the laser spectrum. However, even 
for detunings as large as the width of the pulse spectrum as in 
our example, the phase effect is small and may be neglected 
to first order. Under this approximation, the modulation sig- 
nal can be interpreted as a direct image of the dipole velocity 
after the pump pulse, i.e. the time-evolution of the oscillator 
motion can be retrieved. Further note that our results have 
been obtained for zero CEP of the pulses. While the modu- 
lation parameter Mw(Af) is CEP-independent, the CEP will 
occur as a direct phase offset in the dipole velocity. Hence, 
in the general case the dipole velocity can be reconstructed by 
Vd(0 cos(co r es^ + (pcE)e~ rt up to constant factor, where (Ores 
and 7 are obtained from the LS A fit, and time is measured with 
respect to the pulse peak. 

Coming back to the nonresonant scenario in Fig. [TJl, there 
remains no notable excitation after the pump pulse. The mod- 
ulation parameter Mw(Af) almost exclusively images the triv- 
ial interference-induced fluence variations (interferometric au- 
tocorrelation) of the two pulses at small delays, as can be in- 
ferred from comparison to the M^ CA (At) data, cf. Eq. (13). 



This behavior allows one to identify a nonresonant scenario. 

3 Imaging plasmons in simple-metal clusters 

We now turn to the application of the spectral interferome- 
try scheme to the nonlinear response domain and present a 
theoretical case study on simple-metal clusters. To describe 
the strong- field induced excitation and ionization dynamics 
of the model system Nai47 we apply a semiclassical time- 
dependent density-functional approach on the Vlasov level 
that has b een use d for strong-field laser-cluster interactions 
previousl y^ 3 ! 42 ! 43 !. As test systems we consider the cluster 
ground state and an expanded configuration to study the ef- 
fect of the spectral position of the plasmon resonance. For the 
expanded system we further analyze the dependence of the 
cluster response on pulse intensity. It is demonstrated that the 
spectral interferometry analysis is capable of extracting dy- 
namical observables with sub-fs time resolution and can be 
performed with experimentally accessible observables. 

3.1 Semiclassical Vlasov approach 

The applied Thomas-Fermi- Vlasov approach is a semiclassi- 
cal approximation to time-dependent density-functional the- 
ory and describes the electron dynamics in terms of a contin- 
uous one-particle electron phase-space density /(r,p,f). The 
time evolution of /(r,p,f) is derived from the quantal mean- 
field dynamics by applying the well-known h —> expan- 



which is the equation of motion for the electronic degrees of 
freedom. The effective mean-field potential V G s(r,t) follows 
from variation of the density-dependent total potential energy 
functional with respect to the density as 

Mr) = £Vioa(r-Ri(0) +Vkar(r) + V«(r) +^E(0 -r, 

(17) 

containing the sum over the ion potentials for the present 
configuration R/(f), the electron Hartree potential Vnar, the 
exchange-correlation potential V xc in local-density approxi- 
mation (LDA) [we use the form given in Ref. 48 J, and the inter- 
action with the laser field E(t) in dipole approximation. The 
Hartree term and the exchange-correlation potential are cal- 
culated self-consistently from the actual total electron density 
n e (r,t) = J d 3 p/(r,p,£). To avoid the expensive propagation 
of strongly localized states, only valence electrons are treated 
explicitly. The interaction with nuclei and core electrons is de- 
scribed by a local pseudopotential Vi on ( r ) for sodium 42 . Clas- 
sical dynamics is assumed for the ionic motion. 



sion 



42 44- 47 



. In lowest order this yields the Vlasov equation 

a/ ^' P ' f) = V p /(r,p,0 • V r V eff (r,0 - P • V r /(r,p,f), (16) 
ot m 



Though Eq. ( [16] ) is of classical nature, limited quan- 
tum effects such as exchange and correlation in LDA are 
kept in the effective potential and the initial conditions for 
the distribution function. The latter is determined from 
the Thomas-Fermi ground state according to /°(r,p) = 
2/ (27r/z) 3 0(/? F (r) —p), where 0(jc) is the Heaviside function, 
p F (r) = yj2m\\i — Veff( r )] i s me l° ca l Fermi momentum, and 
jl is the chemical potential to fix the electron number. Propa- 
gation of the initial distribution /°(r,p) by Eq. (16) consti- 
tutes the Thomas-Fermi- Vlasov dynamics. For solving the 
ground- state problem and the time-propagation we apply the 
test particle method in a parallel particle-mesh code with an it- 
erative multigrid Po isson so lver. Further technical details are 
described elsewher e ! 42 ! 47 ! 49 !. 

3.2 Configurations of the model systems 

We consider two different configurations of our test system 
Nai47 . The first one is the ground state with fully relaxed ion 
configuration. This yields an icosahedral structure with highly 
delocalized valence electrons with a mean density close to the 
bulk electron density of sodium (cf. Fig. [2]). The ground state 
has a well localized plasmon resonance at roughly 3 eV and 
will be used as an off-resonant target. 

As a resonant configuration we consider an artificially ex- 
panded cluster, generated by rescaling of the ionic ground 
state geometry. The scaling factor is chosen such that the 
background ion density yields a classical Mie plasmon in res- 
onance with a laser field at 800 nm. Also for this ion config- 
uration, the electronic ground state is calculated (cf. Fig. [2]). 
Note that such expansion can be realized experimentally by 
applying an earlier activation pulse. Though this more realis- 
tic scheme would result in an expanded system that is charged 
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ground state 



expanded 



Fig. 2 Calculated ion and electron distributions of icosahedral 
Nai47 in the ground state and in the expanded configuration (as 
indicated), both in the electronic ground state (see text). The 3D 
images show ion positions as blue spheres and electron isodensity 
surfaces for n e = 4 x 10 21 cm -3 . The expanded configuration 
follows from scaling the ground state ion geometry by a factor of 
1.55, resulting in a classical Mie plasmon energy of about 1.5 eV. 
The ion root-mean- square radii are 8.8 and 13.6 A, respectively. 



and thermally excited, the electronic ground state is used here 
for convenience and clarity. 

3.3 Nonlinear spectral interferometry results 

We now consider laser excitation of the above described clus- 
ter configurations by a pair of intense 3 fs few-cycle pulses 
with zero carrier-envelope phase at 800 nm. As we are aiming 
at the analysis of the pump-pulse-induced collective electron 
motion we begin with an inspection of key observables for ex- 
citation with the pump pulse only, see upper panels of Fig. [3] 

As a nonresonant reference case, pump-only excitation 
of the unexpanded cluster is considered at laser intensity 
/p ea k = 10 13 W/cm 2 , for selected time-dependent observables 
see Fig. [3^. The fully off-resonant excitation yields only 
very low energy absorption and the cluster dipole response 
is a nearly instantaneous polarization that directly follows the 
laser field. The dipole velocity (blue curve) is phase shifted 
by 71 / 2 with respect to the field oscillation of the laser pulse 
(black curve) and nearly vanishes immediately after the pulse. 
The ion and electron distribution of the cluster after the laser 
pulse are basically equal to the initial state (cf. Fig. [3^). 

In addition, two resonant scenarios are depicted in Figs. [3]) 
and [3J:, showing the response of the expanded cluster config- 
urations for pump-only excitation with intensities 7 pea k = 10 12 
and 10 13 W/cm 2 , respectively. In both cases, strong plasmon 
oscillations are excited, as is reflected in the oscillations of the 
dipole velocity after the end of the laser pulse. Note that the 
eigenfrequencies determined from damped oscillator fits are 
very similar (/icOres ~ 1.6eV) for both intensities. 

Most interestingly, for the more intense pump pulse, the 



plasmon oscillations exhibits a substantially longer lifetime of 
?"res = 7.7 fs over the result for the lower intensity scenario of 
2.8 fs, which is a nonlinear effect of the excitation process. 
This behavior can be understood by considering the higher 
cluster ionization achieved with the stronger pulse, see the rel- 
ative ionization (red curves) in Figs.|3]3 and[3]:. For the higher 
intensity, the remaining electron cloud is confined to a vol- 
ume smaller than the ionic frame of the cluster (compare 3D 
plots in Figs. [3}) and[3]:). When oscillating with not-too-large 
amplitude, the electron cloud traverses only the inner cluster 
region where the ion background potential is nearly harmonic, 
which effectively reduces plasmon dephasing by interactions 
with the anharmonic background potential near the cluster sur- 
face. Hence, the nonlinear effect yields an intensity dependent 
increase of the plasmon lifetime, while the plasmon oscilla- 
tion itself remains an effectively linear phenomenon. This is 
analogous to the nonlinear index of refraction (effective linear 
property) of gases or solids. 

Having understood the main physics of the three test cases, 
we now apply the spectral interferometry analysis. Therefore 
the cluster response for pump-probe excitation is calculated 
for systematic delay-scans. For each run the total energy ab- 
sorption and the final ionization is recorded for 200 different 
delays. Modulation signals Mw(At) and Mg(Af) are then cal- 
culated from the absorption and the final ionization as func- 
tion of pulse delay, see Figs. [3]i-[3]f. Before inspecting their 
evolution in detail, we have to motivate the use of these pa- 
rameters for our spectral interferometry analysis, which was 
formulated for a linear absorption cross section. The key idea 
is that the modification of the probe-pulse induced absorption 
and ionization is changed only weakly by the residual coher- 
ent motion of the plasmon, which allows one to linearize the 
absorption and ionization as function of the residual dipole 
velocity amplitude. This approximation will be shown to be 
applicable outside the pulse overlap region. 

Focussing on the nonresonant case (Fig. [3ji) first, both the 
Mw(At) and Mg(Af) signals show oscillations only for very 
short delays within the region of strong pulse overlap. This is 
the typical signature of a nonresonant excitation, in full agree- 
ment with the behavior predicted within the constant cross 



section approximation discussed in Sec. |2.3| The signals thus 
show an image of the interference of the pulses, though signals 
may not be exactly proportional to the field fluence because of 
the nonlinear system response. The nonresonant nature, how- 
ever, can be unambiguously be identified from the modulation 
signals. 

For the resonant scenarios, the modulation parameters show 
strong oscillations well beyond the pulse autocorrelation, see 
Figs. [3^ and |3f. Focussing on the region outside the pulse 
overlap, a nearly perfect mapping of the pump-pulse induced 
dipole velocity into the modulation signals is observed, i.e. 
Mw(At) ~ Mg(Af) oc Vd(0> where v&(t) is the corresponding 



6 | Journal Name, 201 0, [vol],1 -[8] 



This journal is ©The Royal Society of Chemistry [year] 





20 



time [fs] 
10 



20 



0.5 



lpeak=10 13 Wcm-2 

(ground state) 



pulse envelope 



dipole velocity (v d ) 
field strength 




(expanded) 



Lorentz fit of v d : /zw res =1 -58 eV 
t=2.76 fs 
<p=+0.06 7i 





LSAfitof M w : ft<w res =1.60eV 
r=2.61 fs 
<p=+0.05 n 



ISA fit of M Q : /z« res =1 .63 eV 
t=2.62 fs 
<p=-0.23 n 



fit of M Q : fow res =1.59eV 
r=7.67fs 
^=+0.01 K 



10 20 

pulse delay At [fs] 



10 

pulse delay At [fs] 



20 



10 

pulse delay At [fs] 



20 



Fig. 3 Plasmon dynamics and spectral interferometry analysis of Nai47 under intense 3 fs few-cycle pulses at Hcoq = 1.54eV as calculated 
within the semiclassical Vlasov approach (cluster configurations and laser peak intensities as indicated), (a-c) time-evolutions of the laser 
field, dipole velocity (relative to the peak field amplitude) and cluster ionization for pump-only excitation; Lorentz fits of the dipole velocity 
for t > 6fs; 3D snapshot images at t = 6fs. (d-f) systematic delay scans of the modulation signals of the absorption and cluster ionization 
[Mw(Af), Mg(Af)] and corresponding LSA fits for pump-probe excitation. 



dipole velocity for pump-only excitation. In particular, the dif- 
ferent lifetimes are clearly resolved and the damped oscillator 
fits of the dipole signals and the LSA fits of Afw(Af),ATg(Af) 
yield frequencies and lifetimes that agree within a few percent. 
The extracted phase offsets are close to zero, as expected for 
fully resonant excitation. In should be noted that the modu- 
lation signals can be well described by the fit functions, ex- 
cept for short delays in the high intensity scenario, where de- 
viations due to nonlinearities can be found in the delay range 
At = 7. . .9fs, see Fig. [3]f. Otherwise the assumption of a lin- 
earized response seems to be well justified. 

These results support the claim that nonlinear spectral in- 
terferometry with few-cycle pulses opens up an avenue to- 
wards imaging of strong-field induced plasmons in metal- 
clusters with sub-fs time resolution. In addition, as the to- 
be- scanned delay window can be as narrow as a few tens of 
fs, our approach is expected to provide valuable insights into 
the time-evolution of plasmonic properties in expanding clus- 



ters, which would be of high interest for the whole area of 
nanoplasma science. 

4 Conclusions and outlook 

In conclusion, we have explored the applicability of spectral 
interferometry with few-cycle laser pulses as an imaging tool 
for ultrafast collective electron dynamics in the nonlinear re- 
sponse domain. Two simple approximations have been de- 
rived that can be used for fitting the modulation signals of 
energy absorption, ionization, or other (linearizable) observ- 
ables. The performance of the method has been tested with 
a "numerical experiment" on simple-metal clusters, where the 
time-domain analysis of plasmon oscillations with sub-fs reso- 
lution at high modulation contrast has been demonstrated. Our 
results support that both the time-resolved plasmon motion as 
well as key properties such as frequency and lifetime can be 
accurately determined. Moreover, we have shown that strong- 
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field effects, like the increase of the plasmon lifetime due to 
ionization-induced contraction of the electron cloud could be 
resolved. It should be emphasized that the accurate analysis 
of the dynamics has demonstrated by evaluating experimen- 
tally accessible observables, such as total ionization. The un- 
precedented time-resolution of our method promises new fun- 
damental insights into the strong-field induced dynamics in 
finite many particle systems. 
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